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Abstract. We consider the hypothesis that galactic dark matter is composed of ultra- 
i— — i. light scalar particles and use internal properties of dwarf spheroidal galaxies to establish a 

preferred range for the mass of these bosonic particles. We re-investigate the problem of 
the longevity of the cold clump in Ursa Minor and the problem of the rapid orbital decay of 
the globular clusters in Fornax and dwarf ellipticals. Treating the scalar field halo as a rigid 
background gravitational potential and using A-body simulations, we have explored how the 
fS| ' dissolution timescale of the cold clump in Ursa Minor depends on ttia. It is demonstrated that 

■ for masses in the range 0.3 x 10 -22 eV < < 1 x 10~ 22 eV, scalar field dark halos without 

t-h \ self-interaction would have cores large enough to explain the longevity of the cold clump in 

Ursa Minor and the wide distribution of globular clusters in Fornax, but small enough to 
• • ■ make the mass of the dark halos compatible with dynamical limits. It is encouraging to see 

. £h . that this interval of is consistent with that needed to suppress the overproduction of 

substructure in galactic halos and is compatible with the acoustic peaks of cosmic microwave 
radiation. On the other hand, for self-interacting scalar fields with coupling constant A, 
values of m^/X < 0.55 x 10 3 eV 4 are required to account for the properties of the dark halos 
of these dwarf spheroidal galaxies. 
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A Wake created by a gravitational perturber in an infinite homogeneous 



1 Introduction 

The Concordance Cosmological Model, usually referred to as the A+ cold dark matter 
(ACDM) model, has proved to be very successful in explaining observables across a wide rage 
of length scales, from the cosmic microwave background (CMB) anisotropy to the Lyman-a 
forest. In this model, nonbaryonic collisionless cold dark matter (hereafter CDM) makes up 
23% of the total mass content of the Universe. 

Observational data on galactic scales seem to disagree with CDM predictions. High 
resolution A^-body simulations have shown that the predicted number of subhalos is an order 
of magnitude larger than what has been observed [1, 2]. Another discrepancy arises when 
comparing the density profiles of dark halos predicted in simulations with those derived 
from observations of dwarf spheroidal (dSph) galaxies and Low Surface Brightness galaxies 
(LSB's); A r -body simulations predict an universal cuspy density profile, while observations 
indicate that a cored halo is preferred in an important fraction of low-mass galaxies [3-5]. 

These discrepancies might be overcomed by considering other alternatives. One intrigu- 
ing possibility is to consider DM particles as spin zero bosons of ultra- light mass m^, having 
such a large Compton wavelength and a very large number density that they can be described 
by a classical scalar field eft (hereafter Scalar Field Dark Matter; SFDM) [6-20]. In this sce- 
nario the massive scalar field only interacts gravitationally with the rest of the matter and 
is minimally coupled to gravity. Different potentials have been proposed to the scalar field, 



e.g. V($) = V [cosh(£<I») - 1] [13], V($) = m 2 <S> 2 /2 [21, 22] and V($) = m 2 ® 2 /2 + A$ 4 / 4 



[20]. At cosmological scales, it is well known that if the quadratic term of Va, is dominant, 
the SFDM behaves as CDM and thus linear perturbations of SFDM evolve as those in the 



medium 
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standard CDM paradigm [20-22]. However, a difference between SFDM and CDM models is 
that a cut-off in its mass power spectrum, which can prevent the overproduction of substruc- 
tures, naturally arises in SFDM. For a quadratic potential, the substructure overproduction 
issue could be overcome if the the boson associated to the scalar field is ultra-light with a 
mass of ~ 10~ 23 — 10~ 22 [10, 13]. It is important to note that an ultra-light bosonic DM 
particle of mass ~ 1 x 10~ 22 eV is compatible with the acoustic peaks of the CMB radiation 
[23]. 

On the other hand, it has been shown that the scalar field can form stable structures 
that could account for DM halos [24]. It is known [25-31] that ground states are stable 
against spherical and non-spherical perturbations and, furthermore, that such configurations 
are late-time attr actors for quite arbitrary initial profiles of the scalar field [31]. These 
findings suggest that ground states could be naturally formed from initial fluctuations of the 
hypothetical SFDM. 

Interestingly, SFDM halos would have cored mass density profiles. Arbey, Lesgourgues 
& Salati [14] have analyzed the rotation curves in a sample of LSB galaxies with SFDM and 
conclude that ground states can explain fairly well the observed rotation curves, if the mass of 
the boson is ~ 10 -24 — 10 -23 eV (see also [6, 32-34]). If the scalar field is self-interacting, 
its mass and its self-coupling parameter A are both constrained in order to fit the rotation 
curves of spiral galaxies [35-38]. The dynamics of head-on collisions of ground state halos 
have been also studied [39, 40]. 

Here we consider the DM halos of dSph galaxies in the context of SFDM. Our aim is 
to provide new constraints on SFDM assuming that it behaves as a massive and complex 
scalar field. We will discuss cases with and without self-interaction. dSph galaxies provide 
a unique testing ground for the nature of DM. There is growing evidence that Ursa Minor 
(UMi), Fornax, Sculptor, Carina, Leo I and Leo II possess cored DM halos instead of cuspy 
[4, 41-45]. It is worthwhile exploring if the SFDM scenario can explain the existence of 
cores in dSph galaxies and circunvent some problems in galactic dynamics in a natural way. 
Here we investigate two of these problems. First of all, we argue that the interpretation 
that the stellar clump in UMi is a 'dynamical fossil' gives constraints on ms and A. To do 
so, we perform iV-body simulations of the evolution of a cold stellar clump embedded in 
the scalar field DM halo, treating it as a rigid background potential, but including clump's 
self-gravity. Secondly, we argue that, for some combinations of and A, the reduction of 
the gravitational dynamical friction in a SFDM halo could help alleviate the problem of the 
orbital decay of globular clusters (GCs) in dSph galaxies and dwarf ellipticals. 

The article is organized as follows. In §2 we describe the SFDM model and briefly review 
the Schrodinger-Poisson system. In §3, we derive constraints on the values of ms by studying 
the puzzling internal dynamics of UMi dSph galaxy. In §4 we discuss the implications of 
SFDM in the Fornax dSph galaxy. Finally, in section §5 we discuss the results and give our 
conclusions. 



2 The Schrodinger-Poisson System 

Since galactic halos are well described as Newtonian systems, we will work within the New- 
tonian limit. In this limit, the Einstein-Klein-Gordon (EKG) equations for a complex scalar 
field <I> minimally coupled to gravity and endowed with a potential V(3>) = m\<& 2 /2 + A$ 4 /4, 
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can be simplified to the Schrodinger-Poisson equations (SP): 

ihd t iP = -J^V 2 %b + Um^ + -^\^ , (2.1) 

V 2 U = ^Gm\^\ (2.2) 

where is the mass of the boson associated with the scalar field, U is the gravitational 
potential produced by the DM density source (p = 77i^|'0| 2 ), A is the self- interacting coupling 
constant, and the field tft is related to the relativistic field $ through 

[46-48]. 

We are interested in spherical equilibrium solutions to the SP system. They are obtained 
assuming harmonic temporal behavior for the scalar field 

^(r,t) =e-* y */V(r). (2.4) 

In dimensionless variables we have 



m^c m^c 2 



t, (2.5) 



V^Gh U 2 . . 

U -> -z, 7 -> 7770C 7, (2.6) 



c 2 ' c 2 



and the stationary SP system can be written then as 



where 



d 2 (r<ft) = 2r(U - 7) + 2rA^ 3 , (2.7) 
d 2 {rU) = re/) 2 , (2.8) 



A = ^fpA, (2.9) 
87T mi fc 1 



and TTip is the Planck mass. In order to guarantee regular solutions, we require that the 
boundary conditions at r = are d r U = 0, d r (ft = and <^>(0) = <ft c , where (ft c is an arbitrary 
value. On the other hand, we impose that 

POO 

M= <J) 2 r 2 dr, (2.10) 
J 

to be a finite number. 

The system (2.7)-(2.8) can be treated as an eigenvalue problem: given a central value <p c 
and a specific value of A, there are discrete values {7^} for which the solutions {7^, <fti(r), Ui(r)} 
satisfy the boundary conditions. Solutions can be qualitatively differentiated by the number 
of nodes of the radial function (fti . The solution for which (ft has zero nodes is called the ground 
state. The first excited state is the name for the solution for which (ft has one node and so 
on. In this work, we consider only ground states because they are stable under gravitational 
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perturbations and thus are more suitable to model galactic halos than the unstable excited 
states [26]. 

Figure 1 shows the SFDM density profile and the gravitational potential U, as a function 
of radius r, for a ground state with (j) c = 1 and A = 0. In order to specify the properties of a 
certain halo, we will give the mass of the configuration M, and the radius rgs, defined as the 
radius of the sphere that contains 95% of the mass. In particular, the configuration in figure 
1 has a mass M = 2.062 and rgs = 3.799 in dimensionless units (see Eqs. 2.5 and 2.6). 




Figure 1. Density profile (top) and gravitational potential (bottom) of the scalar field in a ground 
state with </> c = 1 and A = 0. 



An interesting characteristic of the SP system is that it obeys the scaling symmetry 

4>, U, 7 , A, r -> eV, e 2 U, e 2 7 , e 2 A, e~V. (2.11) 

This property allows us to find a family of solutions from a particular solution. Because of 
this scaling symmetry, the physical quantities of the solutions satisfy 

M.rgs-^eM.e-W (2.12) 

Therefore, three quantities 4> c , and A or, equivalently e, and A define a model com- 
pletely. It is interesting that for A = 0, cp c ~ 10~ 6 (e ~ 10~ 3 ) and m<* ~ 10 -23 eV, the ground 
state would model a halo with M ~ 1O 1O M0, rgs ~ 7 kpc, and a core radius, r c , defined as 
the radius where the DM density decays a factor of 2 its central value, of ~ 3 kpc. 

While and A are parameters of the SFDM model that once fixed they should be 
universal, <j) c (or equivalently e) is a quantity that may vary from galaxy to galaxy. Indeed, 
it may be directly related to the central density of the DM halo. A typical value of e can be 
inferred from the (dimensionless) compactness of the ground state defined by 

C-*». (2.13) 

^95 

Assuming that dark halos are Newtonian, it is expected that C < 1. This condition is 
fulfilled if e <C 1- In fact, studies of general relativistic equilibrium solutions, show that 
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configurations with a very small central value of the scalar field <3? c have small values of their 
compactness [49]. Furthermore, it is shown that equilibrium configurations with $ c < 10~ 6 
can be treated as Newtonian, that is, they can be described within a good approximation by 
the SP system. Because of the relation (2.3), the condition $ c < 10 -6 implies e < 10 -3 . 

3 The halo of UMi dSph: Constraints on SFDM 

UMi is a diffuse dSph galaxy located at a distance of 69 ± 4 kpc [50] from the Milky Way 
center and has a luminosity of Ly = 3 x 10 5 Lq [50]. Its stellar population is very old with an 
age of 10-12 Gyr. Dynamical studies suggest that UMi is a galaxy dominated by DM, with a 
mass-to- light ratio larger than 60Mq/Lq. Among the most puzzling observed properties of 
UMi is that it hosts a stellar clump, which is believed to be a dynamical fossil that survived 
because the underlying DM gravitational potential is close to harmonic [4, 51]. This condition 
is accomplished if the DM halo has a large core. 

We study whether SFDM halos may have cores large enough to guarantee the survival 
of UMi's clump, but small enough to ensure that the galaxy is not too massive. Indeed, 
models with large cores have large masses. Therefore, we have to impose an upper limit 
on V max , the maximum of the circular velocity of the halo. We will see that this provides 
a stringent constraint on the mass of the boson associated with the scalar field and on the 
self-interaction parameter. 

3.1 The dark halo and the dynamical fossil in UMi 

As we said before, the most remarkable feature in UMi structure is the double off-centered 
density peak [52]. The second peak or clump is located on the north-eastern side of the major 
axis of UMi at a distance of ~ 0.4 kpc from UMi's center. The velocity distribution of the 
stars contained in the clump is well fitted by two Gaussians, one representing the background 
(velocity dispersion of 8.8 km s _1 ) and the other representing the velocity dispersion of the 
second peak (0.5 km s _1 ). The most appealing interpretation is that UMi's clump is a 
disrupted cluster [53] with an orbit in the plane of the sky, which has survived in phase-space 
because the underlying gravitational potential is harmonic [4], implying that the dark halo 
in UMi has a large core. 

In order to explore the parameter space of galactic SFDM halos compatible with the 
survival of the kinematically-cold structure in UMi, we have studied the evolution of a stel- 
lar clump inside a rigid SFDM halo in the ground state. The dynamics of clump's stars 
are simulated with 600 particles in an iV-body code. Hence, a star in the clump feels the 
gravitational force of a static dark halo potential and the gravity of the remainder of the 
stars in the clump. In our simulations we ignore the contribution of the stellar background 
component to the gravitational potential because its mass interior to ~ 0.4 kpc (9 x 10 5 Mq) 
only accounts about 9 percent (or less) the DM mass in UMi. 

We get the gravitational potential U by solving the SP equations. In order to specify a 
SFDM halo we need to fix e and m^. Strigari et al. [54, 55] considered dark halos compatible 
with the observed stellar kinematics of the classical dSph galaxies, including UMi. They 
found that, for realistic density profiles, the mass interior to 300 pc is ~ 1O 7 M for all dSph 
galaxies in the Milky Way halo. Therefore, we fix the e- value in each SFDM halo by imposing 
that the mass within r = 0.39 kpc is 1.5 x 1O 7 M (see also [56-58]), and explore the evolution 
of the clump for halos with different m^. 
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Figure 2. Snapshots of the clump in UMi galaxy, at t = 0, 5, and 10 Gyr. We set the clump on a 
circular orbit in the (x,y)-plane at a distance of r = 0.39 kpc from UMi's center. The mass of the 
boson is = 10~ 23 cV and A = 0. The total mass of the galaxy is M = 9.7 x 1O 9 M . 

Initially, the clump has a density profile p = po exp(— r 2 /2tq) (with r$ ~ 35 pc [59]) 
and was dropped at a galactocentric distance of 0.39 kpc on a circular orbit in the (x, y) 
plane (see also [60]). We explored other eccentric orbits but the survival timescale is very 
insensitive to the eccentricity. If the U-band M/L of UMi stellar population is supposed 
to be Y* = 5.8, the mass of the clump is M c = 7.8 x 1O 4 M . The initial one-dimensional 
velocity dispersion of the clump is 1 km s" 1 [60]. 

The iV-body simulations were performed with the code described in [51] and [60]. This 
iV-body code has the ability to link an individual time step to each particle in the simulation. 
Only for the particle that has the minimum associated time, the equations of motion are 
integrated. This "multi-step" method reduces the typical CPU times of direct particle-particle 
integrations. In our simulations all the particles have the same mass and a softening length 
of 1 pc which is approximately l/10th of the typical separation among clump's particles. The 
convergence of the results was tested by comparing runs with different softening radii and 
number of particles. 

3.2 Scalar field without self- interact ion (A = 0) 
3.2.1 Upper limit on 

Figure 2 shows snapshots of the clump at t = 0, 5 and 10 Gyr for = 10 -23 eV. The orbit 
of the clump lies in the (x, y)-plane. We see that the clump remains intact during one Hubble 
time. This is a consequence that the SFDM halo has a core large enough to guarantee the 
survival of the clump. For a mass of the boson of 10" 23 eV, the core radius is 2.4 kpc and 
the total halo mass of 9.7 x 1O 9 M , which is too massive (see §3.2.2 for a discussion). 

As the shape of the underlying gravitational potential depend on ms, the longevity of 
the clump should depend on m^. Since both the size of the core and the total mass increase 
when ms decreases (see figure 3) , we consider now the evolution of the clump in models with 
nifj, values larger than 10 -23 eV. Figure 4 shows three snapshots of the clump (t = 0, 5, and 
10 Gyr) in a halo with = 10 -22 eV, which has a core radius of ~ 0.7 kpc. The clump 
does not dissolve in one Hubble time. However, when the mass of the boson is 5 x 10~ 22 eV, 
the core radius of the SFDM halo is quite small, of 0.28 kpc, and, consequently, the clump 
loses its identity due to tidal effects (see figure 5). 
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Figure 3. Mass interior to any given radius and density profiles for SFDM halos for models with 
A = and different m^. 
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Figure 4. Same as figure 2 but for a boson mass of = 10 22 eV and A = 0. The UMi's mass is 
of M = 3.1 x 1O 8 M . 

To quantify the destruction of the clump in our simulations, we calculated a map of the 
projected surface density of mass in the (x, y)-plane at any given time t in the simulation. 
We sample this two-dimensional map searching for the 10 x 10 pc size parcel that contains 
the highest mass, This parcel is centered at the remnant of the clump. Figure 6 shows 

the evolution of II with time for models with different and Table 1 summarizes the results 
of the simulations with A = 0. We see that in models with ~ 3 x 10~ 22 eV, the clump 
is diluted within one Hubble time. In halos with > 3 x 10 -22 eV, the clump is erased 
in a too short time. Therefore, we conclude that the survival of the dynamical fossil sets an 
upper limit to the mass of the boson of < 3 x 10 -22 eV. 
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Figure 5. Same as figure 2 but = 5 x 10 22 eV and A = 0. In this case, the total halo mass is 
M = 0.33 x 10 s M n . 
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Figure 6. Temporal evolution of IT for models with A — and different values. 
3.2.2 Lower limit on 

In our SFDM models of UMi's halo, there is a positive correlation between the size of the 
core r c and the maximum of the circular velocity VJnax- While large values of the core are 
favored to explain the persistence of the clump and are allowed by the velocity dispersion 
data, they require very large values of V^ax and M. The total mass of the halos are given in 
Table 1. For instance, with = 10 -23 eV, the total dynamical mass is M = 9.7 x 1O 9 M 
and Vmax — 90 km s" 1 . It is extremely unlikely that one of the less luminous classical dSph 
galaxies has such a large mass [61]. Wilkinson et al. [56] estimate the enclosed mass of UMi 
within ~ 40' (~ 0.8 kpc) to be of the order of 2 x 1O 8 M based on projected radial velocity 
dispersion profiles. Penarrubia et al. [57] suggest that the virial mass of UMi derived from 
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Table 1. Destruction times of the clump in UMi dwarf galaxy for the A = case, for different m^. 
The total mass M, and the radius rgs and r c , are also given for each model. 
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> 10 


1 


11.4 


3.14 


2.13 


0.73 


> 10 


2 


8.3 


1.15 


1.4 


0.50 


> 10 


3 


7.0 


0.65 


1.15 


0.39 


~ 5.0 


4 


6.3 


0.44 


0.96 


0.33 


~ 1.5 


5 


5.9 


0.33 


0.82 


0.28 


~ 0.5 



the assumption of a NFW halo model is ~ 4.5 x 10 9 M Q . In fact, according to the ACDM 
simulations by Zentner & Bullock [62], only 5% of the subhalos in a Milky Way-sized halo has 
total masses > 5 x 1O 9 M0. We will take this value as the upper limit on UMi mass. Note that 
this is a very generous upper limit for M; fits to the velocity dispersion profile give a mass 
of 2.3±} f x 1O 8 M within a tidal radius of 1.5 kpc [54]. If we demand that M < 5 x 1O 9 M , 
we obtain > 1.7 x 10~ 23 eV, whereas for M < 1O 9 M we have > 5 x 10 -23 eV. 
Since there is evidence that at least three dSph galaxies have a dark halo with a large core 
(Fornax, Sculptor and UMi) and ACDM simulations have shown that it is unlikely to have 
a population of three subhalos with M ~ 5 x 1O 9 M in a Milky Way-sized halo (see, e.g., 
[62]), a limit m > 3 x 10" 23 eV seems very reasonable. 

Combining the above lower limit with the upper limit derived in §3.2.1, we find that 
77i0 should lie between 0.3 x 10~ 22 eV to 3 x 10 -22 eV, being our most preferred value 
777^ ~ 1 x 10~ 22 eV. For the latter rriA, values, the core radius of UMi dark halo is between 
500 and 750 pc. 

3.3 Scalar field with self-interaction (A ^ 0) 

So far, we have assumed that boson self- interaction is negligible (A = 0). In order to see how 
7770 depends on self-interaction, we explored models with the third term of Eq. (2.1) being 
distinct from zero (A ^ 0). 

3.3.1 Small A- values (A « 1) 

In this section we consider small values of A (in dimensionless units). The parameters of the 
models are summarized in Table 2. Figures 7 and 8 show the evolution of II(t) for A = 1/2 
and for A = 2, respectively. The halo density profiles are given in figures 9 and 10. With 
A = 1/2, the clump survives even if one increases the mass of the boson up to 4 x 10 -22 eV 



- 9 - 



1 00 




Figure 7. Evolution of II as a function of t for a self-interaction parameter of A = 1/2. 



(here the core radius is ~ 420 pc). In halo models with = 1 x 10~ 21 eV, the clump is 
disrupted in ~ 0.4 Gyr for A = 1/2 but it survives one Hubble time for A = 2. As Table 2 
shows, for the same m^, the core radii of DM halos and their mass increase with A. With 
self-interaction, the permitted window for the mass of the bosonic particles is shifted to 
larger values. We must notice, however, that in order to have M > 10 8 Mq, as derived by 
[54], uiff) < 6 x 10 -22 eV for A = 2. These constraints are compatible with those found in 
[63]. 

3.3.2 Large A-values (A > 1) 

Here, we analyze the effects on the persistence of the clump of UMi for large values of A. It 
was first pointed out in [64] that the inclusion of self-interaction, even for small values of it, 
can lead to important changes in the resulting scalar field configurations. In fact, from the 
definition of the dimensionless variable in equation (2.9), values of cA/h 3 ~ 1 imply A 3> 1. 
In order to study the system of equations (2.7) and (2.8) with A 3> 1, it is useful to rescale 
the variables as follows: (f>* = A 1 / 2 ^ and r* = A _1//2 r. In the limit A — > oo, the Schrddinger 
equation (2.7), to first order in A" 1 , has a simple algebraic form = 7 — U. Substituting 
this relation into the Poisson equation (2.8) and solving it, we obtain an exact solution for 
the scalar density, pj, = 2 , which in terms of the original variables reads as 

_ _ 4>l sin(A-V2 r ) 



A A-V2, 



where <^> 2 is an arbitrary constant to be determined. The condition > provides the 
maximum radius R ma x of the SFDM halo, 



(3.2) 



The mass profile M(r) is given by 
rr 

M(r) = / W = <j) 2 c A 1/2 
Jo 

and the total halo mass M defined as the mass contained within R max reads as 

M = tt^A 1 / 2 . (3.4) 

Prom dimensionless expressions (3.1)-(3.4), it is noticed that the free parameters of the halo 
are A and <p c . However, when the physical dimensions for those expressions are taking into 
account it turns out that the free parameters to fit UMi halo are A 1 / 2 /to 2 and <j> c . We 
constrain these parameters through the following steps 

• assume a reasonable maximum UMi's radius RjjMi- Then, use equation (3.2) to derive 
A 1 / 2 /™ 2 . 

• given the value of A 1 / 2 /m 2 , <p c is obtained from equation (3.3) by requiring that the 
mass (in units of m 2 /m^) contained within a radius of 0.39 kpc is 1.5 x 10 7 Mq. 

• once A 1 / 2 /vrA and 4> c are known, the total halo mass and the core radius are calculated 
and compared with the dynamical limits discussed in §3.2. 

For instance, if we take UMi's tidal radius of 1.5 kpc as the maximum radius of its halo, 
we obtain to^/A ~ 10 3 eV 4 (H = c = 1) and 2 ~ 9.6 x 10 -9 . It follows that the total mass 
is M = 3 x 10 8 Mq and the core radius r c = 900 pc. This mass is consistent with the values 
derived by [54], and the core radius is large enough to guarantee the survival of the clump 
(see tables 1 and 2). We explored different values for the tidal radius of the halo and our 
results are summarized in table 3. We can see that for a tidal radius of ~ 1.8 kpc, the clump 
survives because the core is large enough (~ 1 kpc). In this case, the SFDM halo becomes 
a bit more massive. On the other hand, if we impose the condition M > 1.5 x 10 8 Mq, then 
we obtain to^/A < 1.7 x 10 3 eV 4 . 

4 The halo of Fornax and the orbital decay of GCs 

The fact that SFDM halos have cores might solve other apparent problems in dwarf galaxies. 
In this section we consider the timing problem of the orbit decay of GCs in dwarf elliptical 
galaxies and dSph galaxies. In fact, in a cuspy halo, GCs in these galaxies would have 
suffered a rapid orbital decay to the center due to dynamical friction in one Hubble time, 
forming a nucleated dwarf galaxy. For instance, under the assumption that mass follows 
light or assuming a NFW profile, Fornax GCs 3 and 4, which are at distances to the center 
< 0.6 kpc, should have decayed to the center of Fornax in ~ 0.5-1 Gyr [41, 42] (see also 
Fig. 7 in [65]); this clearly represents a timing problem. Assuming a cuspy NFW halo with 
the same parameters as those reported in [65], GCs 1, 2, 3 and 5 can remain in orbit as 
long as their starting distances from Fornax center are > 1.6 kpc, whereas GC 4 needs an 
initial distance > 1.2 kpc. However, there is no statistical evidence to suggest that the initial 
distribution of GCs is so different to the stellar background distribution. Lotz et al. [66] 
found that the summed distribution of the entire GC systems in 51 Virgo and Fornax cluster 
dwarf ellipticals closely follows the exponential profile of the underlying stellar population. 
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Figure 8. Same as figure 7 but for A = 2.0. 




Figure 9. Mass interior to any given radius and density profiles for SFDM halos for models with 
A = 1/2 and different m^. 

In addition, studies of the radial distribution of GCs in giant elliptical galaxies show that the 
distribution of metal-rich GCs matches the galaxy light distribution [67, 68]. Assuming that 
GCs formed along with the bulk of the field star population in dwarf galaxies, the probability 
that Fornax GCs were formed all beyond 1.2 kpc is ~ (0.03) 5 ~ 2.5 x 1(T 8 (here we have 
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Figure 10. Mass interior to any given radius and density profiles for SFDM halos for models with 
A = 2 and different m^. 



Table 2. Same as Table 1 but for models with small self-interaction parameter. 

A 7710 e M rg5 r c destruction time 

(10" 22 eV) (10~ 5 ) (10 7 M Q ) (kpc) (kpc) (Gyr) 



0.5 



8.13 



17.13 



1.57 0.62 



> 10 



0.5 



5. 



5.53 



4.66 0.92 0.36 



0.5 



10. 



4.68 



2.0 



0.55 0.22 



0.4 



2.0 
2.0 
2.0 
2.0 



2. 

5. 
10. 
20. 



7.9 53.27 1.95 1.06 



5.1 13.73 1.21 0.66 



3.75 5.05 0.82 0.45 



3.0 2.02 0.51 0.2? 



> 10 

> 10 

> 10 
~ 1.0 



used that the fraction of the stellar mass beyond 1.2 kpc for a King model with a core radius 
of 0.4 kpc and a concentration parameter of 0.72, as used in [65], is less than 3 percent). 
Therefore, it is very unlikely that all the GCs in Fornax were formed at such large distances 
and even they did, there is still a timing problem with GCs 3 and 4. 
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Figure 11. Left panel: Circular velocity (solid line) and sound speed (dotted line) as a function of 
galactocentric distance for a SFDM model with ml /A = 0.55 x 10 3 eV 4 and i? max = 2.1 kpc. For 
comparison we also plot the circular velocity in the cored collisionless DM halo employed by Sanchez- 
Salccdo ct al. [42] to solve the problem of orbit surviving. Right panel: Gravitational drag force in 
these models. Solid line stands for the SFDM model and dashed line for the the cored collisionless 
halo. 



The dynamical friction timescale depends on the nature of DM particles. Goodman [69] 
was the first to suggest that the suppression of dynamical friction in superfluid DM halos 
could circumvent some problems in galactic dynamics such as the presence of rotating bars in 
dense dark halos (see also [70, 71]). Here we will study if SFDM halos can solve or alleviate 
the timing problem of the GCs in dwarf galaxies. 

We consider first the so-called Thomas- Fermi regime, that is, the limit of large scattering 
parameter (A ^> 1). In this limit, the dynamics of the scalar field is effectively the same as 
that for an ideal fluid with equation of state P = Kp 2 and K = A/m^ (see [71] and the 
Appendix). Therefore, the SFDM dynamical friction force exerted on a gravitational object 
is the same as it is in a gaseous medium with sound speed c s = (2Kp) 1 / 2 (see the Appendix). 

Ostriker [72] calculated the drag force on a perturbing body of mass M p moving at 
constant speed Vq through an ideal fluid with sound speed c s . She found that for subsonic 
perturbers, that is Ai < 1 where the Mach number is defined as M = Vo/c s , the drag force 
is generally smaller in a gaseous medium than in a collisionless medium, because pressure 
forces in a gaseous medium create a symmetric distribution in the vicinity of the perturber. 
More specifically, she found that at M < 1 the dynamical friction force Fdf hi a gaseous 
and homogeneous medium is 



F DF 



Att{GM p Y Po 



v 2 



2 \1-M 



(4.1) 



where po is the unperturbed fluid density. Moreover, Kim & Kim [73] demonstrated that 
Ostriker's formula for the subsonic case, which was derived for perturbers in rectilinear orbits, 
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is also valid for objects on circular orbits. 

We consider now the decay of GCs on near-circular orbits in a spherical self-gravitating 
halo of fluid DM with polytropic index 7 = 2, due to dynamical friction. In that halo, the 
circular speed depends on the galactocentric distance r, Vo(r) = y / GM(r)/r, where M(r) is 
given in Eq. (3.3). The motion is subsonic, i.e. Vq/c s < 1, at any radius interior to 0.84r c . 
Figure 11 shows the halo circular velocity, sound speed and the drag force (using Eq. 4.1) as 
a function of the distance to the center of Fornax for a model with i? max = 2.1 kpc, which 
corresponds to r c = 1.25 kpc and ra^jX = 0.55 x 10 3 eV 4 . For comparison, we also plot the 
corresponding curves for a King halo made of collisionless particles with a core large enough 
to solve the timing problem in Fornax (see [42] for details). Within a galactocentric distance 
of 0.8 kpc, the dynamical friction in a collisionless halo with a King radius of 1.5 kpc is 
similar to the drag force in a SFDM halo with a core radius of 1.25 kpc and thus, conclude 
that models with va^jX < 0.55 x 10 3 eV 4 are compatible with the radial distribution of the 
population of GCs in Fornax. 

In the noninteracting case, the equation governing the evolution of the perturbation is 
no longer a wave equation with constant sound speed but, instead, the velocity group depends 
on the wavenumber [74] . The derivation of the steady density wake in a homogeneous medium 
is outlined in the Appendix. For a perturber moving at speed Vq in rectilinear orbit along 
the z-axis, it is convenient to define the transonic wavenumber as &o = m^Vo/h and the 
radial wavenumber = k\ + ky. In the Appendix, we show that Fourier modes with 
radial wavenumber kn larger than k$ behaves subsonically, whereas modes with kn < ko are 
supersonic. Consider a GC on circular orbit with radius -R rb = 0.8 kpc from the galactic 
center of Fornax with a circular velocity ~ 18 km s _1 (see Fig. 11). For = (0.1 — 1) x 10~ 22 
eV, the transonic wavenumber is k$ = 0.1 — 1 kpc -1 . Small ko values imply that the motion 
is supersonic only for Fourier modes with large wavelength. Since GCs are assumed to be 
on quasicircular orbits, the direction of Vq changes by tt/2 four times per orbit and thus 
the forward- wave propagation effectively restarts as well [72, 73, 75]. This implies that 
modes with wavelength much larger than R OT h are not relevant. In fact, because the wake 
bends behind the body, the drag force promptly saturates to a steady-state value within 
less than the crossing time of modes over the distance equal to the orbital diameter [73]. 
Therefore, we may state that the perturber is essentially subsonic provided that k^ 1 > i? or b- 
As dynamical friction is effectively suppressed for subsonic perturbers, Fornax GCs within 
a radius of ~ 1 kpc are expected to avoid significant orbital decay if ko < 1 kpc" 1 , which 
implies < 1 x 10~ 22 eV. Thus, values between 0.3 x 10~ 22 eV and 1 x 10~ 22 eV could 
explain the survival of the dynamical fossil in UMi and may alleviate the timing problem of 
the GCs orbiting dwarf galaxies. 

5 Concluding Remarks 

We have considered a model where ultra-light bosons are the main components of the dark 
halos of galaxies. The main goal of this work was to constrain the mass of the scalar particles. 
We constructed stable equilibrium configurations of SFDM in the Newtonian limit to model 
the DM halo in UMi. We studied two relevant cases of SFDM halos: with and without 
self-interaction. 

The persistence of cold substructures in UMi places upper limits on m^. Using N- 
body simulations, we found that the survival of cold substructures in UMi is only possible if 
m^<3x 10 -22 eV in the A = case. On the other hand, by imposing a plausible upper 
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Table 3. Relevant parameters of the halo for models with large self- interaction. 



(kpc) 


(10 3 eF 4 ) 


9 c 

(io- 9 ) 


(kpc) 


l\l 

(10 8 M ) 


1.2 


1.7 


6.35 


0.72 


1.6 


1.3 


1.4 


7.35 


0.78 


2.0 


1.4 


1.2 


8.4 


0.84 


2.5 


1.5 


1.1 


9.6 


0.90 


3.0 


1.6 


0.94 


10.8 


0.96 


3.6 


1.7 


0.83 


12.1 


1.02 


4.3 


1.8 


0.74 


13.5 


1.08 


5.1 



limit on M, we place lower limits on m^. All together, we find that for A = 0, ms should be 
in the window 

0.3 x 10~ 22 eV < < 3 x 10~ 22 eV. (5.1) 

Since the timing problem of the orbital decay of the GCs in Fornax can be alleviated if 
m^<lx 10~ 22 eV for A = 0, our most favored value is around (0.3 — 1) x 10 22 eV. 

For SFDM models with self-interaction, the upper limit on increases with A. For 
A = 2, halos made up by bosons of mass < 6 x 10~ 22 eV could account for the observed 
internal dynamics of UMi. In the limit A > 1, we find that m^/X < 0.55 x 10 3 eV 4 would 
explain the longevity of UMi's clump and the surviving problem of GCs in Fornax. 

The window of permitted values for is quite narrow. Even so, it is remarkable that 
our preferred range for the mass of the boson derived from the dynamics of dSph galaxies is 
compatible with those given by other authors to ameliorate the problem of overabundance 
of substructure and is consistent with the CMB radiation [10, 13, 23]. 

In a recent posting during the course of submitting this paper, Slepian and Goodman 
[71] constrain the mass of bosonic DM using rotation curves of galaxies and Bullet Cluster 
measurements of the scattering cross section of self-interacting DM (non bosonic) under 
the assumption that these systems are in thermodynamic equilibrium. If their assumptions 
are confirmed to be valid, repulsive bosonic DM will be excluded and, thereby, the only 
remaining window open is non-interacting bosons with masses in the range given in equation 
(5.1). Nevertheless, the static diffusive equilibrium between Bose-Einstein condensate and 
its non-condensated envelope, as well as finite temperature effects need to be reconsidered 
[76]. In addition, other authors [63, 77] argue that scattering cross sections for bosonic DM 
are much smaller than those derived from the condition of thermodynamic equilibrium by 
Slepian and Goodman [71]. 
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A Wake created by a gravitational perturber in an infinite homogeneous 
medium 

In order to derive the dynamical friction felt by a gravitational perturber moving through 
a SFDM medium or Bose-Einstein condensate, it is convenient to use the hydrodynamical 
representation [74, 78]. In the quantum-mechanical equation of motion (2.1), we write the 
wave function in the form ip(r, t) = A(r, t) exp(iS(r, t) /K) where A and S are real functions, 
and make the Madelung transformation p = m^\ip\ 2 = m^A 2 and v = ViS/m^,, to find 

^+V(p-v) = 0, (A.l) 

d v 

— + ( v ■ V)v = -V (U + U cxt + Q) - -VP S , (A.2) 
at p 

where P s is the pressure arising from the short-range interactions, U ex t is the moving and 
rigid gravitational potential created by the perturber (a GC in our case) and 

Q = -^?*^ (A - 3) 

is the quantum potential, which arises from the Heisenberg uncertainty principle. The equa- 
tion of state is polytropic, P s = Kp J , with index 7 = 2 and K = A/(4m^) [69, 74]. Thus, we 
can define the sound speed as c s = {2Kp) 1 ^ 2 . 

The Thomas-Fermi approximation amounts to neglecting the quantum potential in 
Eq. (A.2). In that case, Equations (A.1)-(A.2) reduce to the usual Euler equations of a 
fluid. Therefore, all we know about dynamical friction in gaseous media holds valid for a 
Bose-Einstein condensate. 

In the noninteracting case, the particles interact only via gravity and thus P s = 0. We 
study the completely steady flow created by a gravitational perturber moving in a straight-line 
trajectory at constant-speed Vq through a homogeneous (infinite) medium. The gravitational 
attraction between the perturber and its own wake is the dynamical friction force. Consider 
a particle at the origin of our coordinate system, surrounded by a scalar field medium whose 
velocity far from the particle is Vq = —Vqz, with Vo > 0. We will assume that the perturber 
induces a small perturbation and hence describe the medium's response in linear theory. 
Once linearized, the quantum Euler equations (A.1)-(A.2) with P s = are 

p V • v' + V • Vp' = 0, (A.4) 

h 2 



V ■ W = -Wext + ^2 V 



V 2 l 

Po 



(A.5) 
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where p' = p—po and v' = v — Vq. By substituting Eq. (A. 4) into the divergence of Eq. (A. 5), 
we find a single differential equation for the density perturbation, a = p'/po, 

- V ^ = - V2U ^ + ^ V2 ^ a ^ (A - 6) 

which is a three dimensional version of the Bernoulli-Euler beam equation. For a point-mass 
M p , we have that p eyi t = M p 8(x)S(y)6(z) and, hence, V 2 C/ ext = 4irGM p 5(x)5(y)5(z). Taking 
the Fourier transform of Equation (A. 6), we have 

(/c 4 - Aklkl) a = JIgM p , (A.7) 

<t> 

where k$ = m^Vo/h. The density perturbation is recovered by doing the inverse Fourier 
transform, 

2GM v m\ f exp(ifc • r) 3 
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In order to find the poles of the integrand, it is convenient to write the above equation as 



_ 2GMgi , exp(|_r) 



where x± = ± y ~~ ^ anc ^ k 2 R = k 2 + k 2 . The integral (A.9) along & 2 is evaluated by 
transforming to the complex plane. For kn > ko, all the poles are outside the contour of 
integration, whereas the integrand has poles on the real axis for ku < ko- This reflects the 
fact that Fourier modes with kn > ko has a velocity group larger than the perturber's velocity 
and thus are subsonic, whereas modes with small wavenumber (kji < ko) are supersonic 1 . It 
is convenient to subdivide the interval of integration over kn into two parts: [0, ko] (interval 
I) and [A;o,oo] (interval II). Thus a(r) = ai(r) + a n (r), where «i(r) corresponds to the 
fcij-integral over interval I and a u (r) is the integral over interval II. The component a n (r) 
displays front-back symmetry, whereas the contribution cti(r) is nonzero only behind the body 
in order to preserve causality. In terms of the dimensionless distance along the symmetry 
axis z = koz, and the dimensionless cylindrical radius it! = koR, the density perturbation can 
be written as 

a(S, R) = k& R) + Us, R)} , (A.10) 

where £i(z, R) and £11(2, R) are dimensionless quantities, and the subscripts I and II refer to 
the asymmetric and the symmetric parts, respectively. We see that the density enhancement 
at (z,R) decreases with increasing Vo and increases with ko- 

For illustration, the symmetric part of the density perturbation, £ n , is shown in Figure 
12. This map represents the symmetric contribution of the steady-state density perturbation 
created by a gravitational perturber moving in an infinite homogeneous medium. We see 
that the density isocontours are extremely spherical (for a comparison with the isocontours 
in a gaseous medium see [72]). The time-dependent evolution of the wake generated by, and 
the gravitational drag force on, a body traveling on a circular orbit will be presented in a 
forthcoming paper. 



References 



[1] Klypin, A., Kravtsov, A., Valenzuela, O. & Prada, F., Where Are the Missing Galactic 
Satellites?, Astrophys. J. 522 (1999) 82 [astro-ph/9901240]. 

[2] Ostriker, J. P., & Steinhardt, P., New Light on Dark Matter, Science 300 (2003) 1909 
[astro-ph/0306402]. 

[3] van den Bosch, F. C, Robertson, B. E., Dalcanton, J. J. & de Blok, W. J. G., Constraints on 
the Structure of Dark Matter Halos from the Rotation Curves of Low Surface Brightness 
Galaxies, Astron. J. 119 (2000) 1579 [astro-ph/9911372]. 

[4] Kleyna, J. T., Wilkinson, M. I., Gilmore, G. & Evans, N. W., A Dynamical Fossil in the Ursa 
Minor Dwarf Spheroidal Galaxy, Astrophys. J. 588 (2003) L21 [astro-ph/0304093]. 

[5] de Blok, W., J., G. & Bosma, A., High-resolution rotation curves of low surface brightness 
galaxies, Astron. and Astrophys. 385 (2002) 816 [astro-ph/0201276]. 

[6] Sin S.J., Late-time phase transition and the galactic halo as a Bose liquid, Phys. Rev. D 50 
(1994) 3650 [hep-ph/9205208]. 

See [74] for the derivation of the velocity group. 



- 19 - 



[7] Peebles, P.J.E., & Vilenkin, A., Noninteracting dark matter, Phys. Rev. D. 103506 (1999) 60 
[astro-ph/9904396]. 

[8] Peebles, P.J.E., Fluid dark matter, Astrophys. J. L127 (1999) 543 [astro-ph/9904396]. 

[9] Salmi, V. and Wang, L., New cosmological model of quintessence and dark matter, Phys. Rev. 
D 62 (2000) 103517 [astro-ph/9910097]. 

[10] Hu W., Barkana, R., & Gruzinov A., Fuzzy Cold Dark Matter: The Wave Properties of 
Ultralight Particles, Phys. Rev. Lett. 85 (2000) 1158 [astro-ph/0003365]. 

[11] Matos, T., Guzman F. S. & Ureha-Lopez, L. A., Scalar field as dark matter in the universe, 
Class. Quant. Grav. 17 (2000) 1707 [astro-ph/9908152]. 

[12] Matos, Urena-Lopez, L. A., LETTER TO THE EDITOR: Quintessence and scalar dark matter 
in the Universe, Class. Quant. Grav. 17 (2000) L75 [astro-ph/0004332]. 

[13] Matos, T., & Urena-Lopez, L. A., Further analysis of a cosmological model with quintessence 
and scalar dark matter, Phys. Rev. D 63 (2001) 063506 [astro-ph/0006024]. 

[14] Arbcy A., Lcsgourgucs, J., & Salati, P., Quintessential halos around galaxies, Phys. Rev. D 64 
(2001) 123528 [astro-ph/0105564]. 

[15] Lee, Jae-Weon, Is Dark Matter a EEC or Scalar Field?, J. Korean Phys.Soc. 54 (2009) 2622 
[arXiv:0801.1442]. 

[16] Marsh, D.J.E., & Ferreira, P. G., Ultra-Light Scalar Fields and the Growth of Structure in the 
Universe, Phys. Rev. D. 82 (2010) 103528 [astro-ph/0812.1342v3]. 

[17] Urena-Lopez, L. A., Bose-Einstein condensation of relativistic Scalar Field Dark Matter, J. 
Cosmol. Astropart. Phys. JCAP01 (2009)014 [gr-qc/0806.3093]. 

[18] Lundgren, A. P., Bondarescu, M., Bondarescu, R.,& Balakrishna, J., Lukewarm dark matter: 
Bose condensation of ultralight particles, Astrophys. J. 715 (2010) L35 [astro-ph/1001.0051v2]. 

[19] Hwang, J. C, Roles of a coherent scalar field on the evolution of cosmic structures, Phys. Lett. 
B 401 (1997) 241 [astro-ph/9610042]. 

[20] Matos, T. & Suarez, A., Structure formation with scalar-field dark matter: the fluid approach, 
Mon. Not. Roy. Astron. Soc. 416 (2011) 87 [arXiv: 1101. 4039]. 

[21] Turner, M. S. Coherent scalar-field oscillations in an expanding universe, Phys. Rev. D 28 
(1983) 1243. 

[22] Matos, T., Vazquez- Gonzalez, A., & Magaha, J., </> 2 as dark matter, Mon. Not. Roy. Astron. 
Soc. 393 (2009) 1359 [astro-ph/0806.0683]. 

[23] Rodrfgucz-Montoya, I., Magaha, J., Matos, T., & Perez-Lorenzana, A., Ultra Light Bosonic 
Dark Matter and Cosmic Microwave Background, Astrophys. J. 721 (2010) 1509 
[arXiv:0908.0054]. 

[24] Chavanis, P.H. Mass-radius relation of Newtonian self- gravitating Bose-Einstein condensates 
with short-range interactions. II. Numerical results, Phys. Rev. D 84 (2011) 043532, 
[arXiv:1103. 1103.2054]. 

[25] Gleiser, M., Stability of boson stars, Phys. Rev. D. 38 (1988) 2376. 

[26] Lee, T.D., & Yang Pang, Stability of mini-boson stars, Nuclear Physics B 315 (1988) 477. 

[27] Kusmartsev, F.V., Mielke, E.W., & Schunk, F.E., Stability of neutron and boson stars: a new 
approach based on catastrophe theory, Phys. Lett. A 157 (1991) 465. 

[28] Seidcl, E., & Suen Wai-Mo, Dynamical evolution of boson stars: perturbing the ground state, 
Phys. Rev. D. 42 (1990) 384. 

[29] Jctzer, P., Boson stars, Phys. Reports 220 (1992) 163. 



- 20 - 



Balakrishna, J., Bondarescu, R., Daues, G., Guzman F.S., Seidel, E., Evolution of 3D Boson 
Stars with Waveform Extraction, Class. Quant. Grav. 23 (2006) 2631. 

Bcrnal, A., & Guzman, F. S., Scalar field dark matter: Nonspherical collapse and late-time 
behavior, Phys. Rev. D 74 (2006) 063504 [astro-ph/0608523]. 

Ji S. U., & Sin Sang- Jin, Late-time phase transition and the galactic halo as a Bose liquid. II. 
The effect of visible matter, Phys. Rev. D 50 (1994) 3655 [arXiv:hep-ph/9409267]. 

Bernal, A., Matos T., & Nunez, D., Flat Central Density Profiles from Scalar Field Dark 
Matter Halos, Rev. Mex. Astron. Astrof. 44 (2008) 149 [astro-ph/0303455]. 

Jae-Weon Lee, Sz Sooil Lim, Minimum mass of galaxies from BEC or scalar field dark matter, 
JCAP 1001 (2010) 007 [astro-ph/0812.1342v3]. 

Jae-weon Lee <fc In-guy Koh, Galactic Halos As Boson Stars, Phys. Rev. D. 53 (1996) 2236 
[hep-ph/9507385v2]. 

Arbcy A., Lesgourgucs, J., & Salati, P., Galactic halos of fluid dark matter, Phys. Rev. D 68 
(2003) 023511 [astro-ph/0301533]. 

Bohmcr, C. G., & Harko, T., Can dark matter be a Bose Einstein condensate?, J. Cosmol. 
Astropart. Phys. JCAP06(2007)025 [astro-ph/0705.4158v4]. 

Harko, T., Bose-Einstein condensation of dark matter solves the core/cusp problem, J. Cosmol. 
Astropart. Phys. JCAP05(2011)022 [arXiv:1105.2996]. 

Bcrnal, A., &: Guzman, F. S., Scalar field dark matter: Head-on interaction between two 
structures, Phys. Rev. D 74 (2006) 103002 [astro-ph/0610682]. 

Gonzalez, J. A. & Guzman, F. S., Interference pattern in the collision of structures in the 
Bose-Einstein condensate dark matter model: Comparison with fluids, Phys. Rev. D 83 (2011) 
103513 [arXiv:1105.2066]. 

Goerdt, T, Moore, B., Read, J. I., Stadel, J. <fc Zemp, M., Does the Fornax dwarf spheroidal 
have a central cusp or core?, Mon. Not. Roy. Astron. Soc, 368 (2006) 1073 [astro-ph/0601404]. 

Sanchcz-Salcedo, F. J., Rcycs-Iturbide, J. & Hernandez, X., An extensive study of dynamical 
friction in dwarf galaxies: the role of stars, dark matter, halo profiles and MOND, Mon. Not. 
Roy. Astron. Soc. 370 (2006) 1829 [astro-ph/0601490]. 

Battaglia, G., Hclmi, A., Tolstoy, E., Irwin, M., Hill, V. & Jablonka, P., The Kinematic Status 
and Mass Content of the Sculptor Dwarf Spheroidal Galaxy, Astrophys. J. 681 (2008) L13 
[arXiv:0802.4220]. 

Amorisco, N. C. & Evans, N. W., Dark Matter Cores and Cusps: The Case of Multiple Stellar 
Populations in Dwarf Spheroidals, 2011, [arXiv:1106.1062vl]. 

Gilmore, G., Wilkinson, M. I., Wyse, R. F. G., Kleyna, J. T., Koch, A., Evans, N. W & 
Grcbcl, E. K., The Observed Properties of Dark Matter on Small Spatial Scales, Astrophys. J. 
663 (2007) 948 [astro-ph/0703308vl] 

Friedberg, R., Lee, T.D., & Pang, Y., Mini-soliton stars, Phys. Rev. D 35 (1987) 3640. 

Chi-Wai Lai PhD Thesis. http://laplace.phas.ubc.ca/Members/matt/Doc/Theses/Phd/lai.pdf 

A. Bcrnal PhD Thesis, http://www.fis.cinvestav.mx/thesis/indcx.php/thesis 

Guzman, F. Siddhartha & Urcha-Lopcz, L. Arturo, Evolution of the Schrddinger-N ewton 
system for a self- gravitating scalar field, Phys. Rev. D 69 (2004) 124033 [gr-qc/0404014]. 

Grebel, E. K., Gallagher, J. S., Harbeck, D., The Progenitors of Dwarf Spheroidal Galaxies, 
Astronomical. J 125 (2003) 1926 

Lora, V., Sanchez-Salcedo, F. J., Raga, A. C. & Esquivel, A., An Upper Limit on the Mass of 



-21- 



the Black Hole in Ursa Minor Dwarf Galaxy, Astrophys. J. 699 (2009) L113 [arXiv:0906.0951]. 

[52] Kleyna, J. T., Gcller, M. J., Kenyon, S. J., Kutz, M. J., Thorstcnscn, J. R., A V and I CCD 
Mosaic Survey of the Ursa Minor Dwarf Spheroidal Galaxy, Astrophys. J. 115 (1998) 2359. 

[53] Read, J. I., Wilkinson, M. I., Evans, N. W., Gilmorc, G., Kleyna, J. T., The importance of 
tides for the Local Group dwarf spheroidals, Mon. Not. Roy. Astron. Soc. 367 (2006) 387 
[astro-ph/0511759]. 

[54] Strigari, L. E., Bullock, J. S., Kaplinghat, M., Diemand, J., Kuhlen, K & Madau, P., 

Redefining the Missing Satellites Problem, Astrophys. J. 669 (2007) 676 [arXiv:0704.1817]. 

[55] Strigari, L. E., Bullock, J. S., Kaplinghat, M., Simon, J. D., Geha, M., Willman, B., Walker, 
M. G., A common mass scale for satellite galaxies of the Milky Way, Nature 454 (2008) 1096 
[arXiv:0808.3772]. 

[56] Wilkinson, M. I., Kleyna, J. T., Evans, N. W., Gilmore, G. P., Irwin, M. J., Grebel, E. K., 

Kinematically Cold Populations at Large Radii in the Draco and Ursa Minor Dwarf Spheroidal 
Galaxies, Astrophys. J. 611 (2004) 21 [arXiv:astro-ph/0406520]. 

[57] Penarrubia, J., McConnachie, A. W., Navarro, J. F., The Cold Dark Matter Halos of Local 
Group Dwarf Spheroidals, Astrophys. J. 672 (2008) 90 [arXiv:astro-ph/0701780]. 

[58] Walker, M. G., Mateo, M., Olszewski, E. W., Penarrubia, J., Wyn E. N., Gilmorc, G., A 
Universal Mass Profile for Dwarf Spheroidal Galaxies?, Astrophys. J. 704 (2009) 1274 
[arXiv:0906.0341]. 

[59] Palma, C, Majewski, S. R., Siegel, M. H., Patterson, R. J., Osthcimer, J. C. & Link, R., 

Exploring Halo Substructure with Giant Stars. IV. The Extended Structure of the Ursa Minor 
Dwarf Spheroidal Galaxy, Astron. J. 125 (2003) 1352 [astro-ph/0205194]. 

[60] Sanchcz-Salcedo, F. J., & Lora, V., Bounds on the Mass and Abundance of Dark Compact 
Objects and Black Holes in Dwarf Spheroidal Galaxy Halos, Astrophys. J. 658 (2007) L83 
[astro-ph/0702502]. 

[61] Strigari, L. E., Bullock, J. S., Kaplinghat, M., Kravtson, A., V., Gncdin, O. Y., A Large Dark 
Matter Core in the Fornax Dwarf Spheroidal Galaxy?, Astrophys. J. 652 (2006) 306 
[astro-ph/0603775]. Strigari, L. E., Bullock, J. S., Kaplinghat, M., Kravtson, A., V., Gncdin, 
O. Y., A Large Dark Matter Core in the Fornax Dwarf Spheroidal Galaxy?, Astrophys. J. 652 
(2006) 306 [astro-ph/0603775]. 

[62] Zcntner, A. R. & Bullock, J. S., Halo Substructure and the Power Spectrum, Astrophys. J. 598 
(2003) 49 [astro-ph/0304292]. 

[63] Rindlcr-Dallcr, T., Shapiro, P. Angular Momentum and Vortex Formation in 

Bose- Einstein- Condensed Cold Dark Matter Haloes, submitted to MNRAS (2011) 
[arXiv:1106.1256]. 

[64] Colpi, M., Shapiro, S. L., Wasscrman, I., Boson stars - Gravitational equilibria of 
self-interacting scalar fields, Phys. Rev. Lett. 57 (1986) 2485. 

[65] Angus, G. W., Diaferio, A. Resolving the timing problem of the globular clusters orbiting the 
Fornax dwarf galaxy, Mon. Roy. Astron. Soc. 396 (2010) 887. 

[66] Lotz, J. M., Telford, R., Ferguson, H. C, Miller, B. W., Stiavelli, M., Mack, J. Dynamical 
Friction in DE Globular Cluster Systems, Astrophs. J. 552 (2001) 572 

[67] Harris, W. E., Pctrie, P. L. The space distribution of globular clusters in M4-9, Astrophys. J. 
223 (1978) 88. 

[68] Strader, J., Romanowsky, A., Brodic, J., Spitlcr, L., Beasley, M., Arnold, J., Tamura, N., 
Sharpies, R., Arimoto, N. Wide-Field Precision Kinematics of the M87 Globular Cluster 
System, (2004), [arXiv:1110.2778]. 



-22- 



[69] Goodman, J. Repulsive dark matter, New Astron. 5 (2000) 103 [arXiv:astro-ph/0003018]. 

[70] Goodman, J. & Slcpian, Z. Chance and Chandra, Pramana 77 (2011) 107 [arXiv: 1104.2620]. 

[71] Slepian, Z., & Goodman, J. Ruling Out Bosonic Repulsive Dark Matter, arXiv:1109.3844 

[72] Ostriker, E. Dynamical Friction in a Gaseous Medium, Astrophys. J. 513 (1999) 252 
[arXiv:astro-ph/9810324]. 

[73] Kim, H. & Kim, W. T. Dynamical Friction of a Circular- Orbit Perturber in a Gaseous 
Medium, Astrophys. J. 665 (2007) 432 [arXiv:0705.0084]. 

[74] Chavanis, P.H. Mass-radius relation of Newtonian self- gravitating Bose-Einstein condensates 
with short-range interactions. I. Analytical results, Phys. Rev. D 84 (2011) 043531, 
[arXiv:l 103.2050]. 

[75] Sanchcz-Salcedo, F. J., Brandenburg, A. Dynamical friction of bodies orbiting in a gaseous 
sphere, Mon. Roy. Astron. Soc. 322 (2001) 67. 

[76] Harko, T., Madarassy, E. J. M. Finite temperature effects in Bose-Einstein Condensed dark 
matter halos, (2011), [arXiv:1110.2829]. 

[77] Rodn'guez-Montoya I., Perez-Lorenzana, A., De La Cruz-Burelo, E., Giraud-Heraud, Y., 
Matos, T. Cosmic Bose dark matter,(2011) [arXiv:1110.2751]. 

[78] Pitaevskii & Stringari Bose Einstein Condensation,(2004) , Oxford University Press. 



- 23 - 



